Electronic Control of Spin Alignment in vr-Conjugated Molecular Magnets 
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Intramolecular spin alignment in 7r-conjugated molecules is studied theoretically in a model of 
a Peierls-Hubbard chain coupled with two localized spins. By means of the exact diagonalization 
technique, we demonstrate that a spin singlet [S = 0) to quartet [S = 3/2) transition can be 
induced by electronic doping, depending on the chain length, the positions of the localized spins, 
and the sign of the electron-spin coupling. The calculated results provides a theoretical basis for 
understanding the mechanism of spin alignment recently observed in a diradical donor molecule. 

PACS numbers: 75.50.Xx, 71.10.Fd, 75.20.Hr 
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Molecular magnetism has captured a large number of 
experimental and theoretical investigations for decades 
Qi M IE ^ IE IS ' Recently Izuoka and his coworkers have 
succeeded in controlling intramolecular spin alignment by 
charge doping in a newly designed organic molecule, thi- 
anthrene bis(nitronyl nitroxide) 0. This molecule, con- 
sisting of TT-conjugated moiety and two stable radicals, 
is spin singlet [S = 0) in its ground state and becomes 
spin quartet [S = 3/2) upon one-electron oxidation, i.e. 
one-hole doping. The control of spin alignment by pho- 
toexcitation has also been demonstrated in a similar type 
of molecules . The observation of such controllable in- 
tramolecular spin alignment by external stimuli offers us 
very alluring prospect for its future application on ultra- 
fast light-switching or magnetic storage, because of flex- 
ibility and variety in designing 7r-conjugated molecules. 

The dominant mechanism of spin alignment in the 
ground state of magnetic molecules has been well estab- 
lished as a topological rule 0, 0, Q| , based on the dy- 
namical antifcrromagnetic spin polarization effect of tt 
electrons with on-site Coulomb repulsion. However, the 
topological rule cannot be applied to the doped or excited 
molecules, and most discussions on the spin alignment 
control so far have been give n on the basis of individual 
molecular orbitals H EUll H IH . The purpose of the 
present paper, therefore, is to provide a theoretical basis 
for the spin alignment control by charge doping from a 
more general point of view. The central question is how 
the one-electron addition to or removal from the inter- 
acting TT-electron system alters the correlation between 
the two separate localized spins. 

For this purpose, we design a theoretical model to in- 
corporate delocalized 7r-electrons and localized spins: a 
Peierls-Hubbard model on a A'^-site linear chain with two 
localized spins interacting with the tt electrons in the 
form of the Kondo coupling. The linear chain part cor- 
responding to polyene is subject to electron-lattice in- 
teractions, which are taken into account in the form of 
the Su-Schrieff-Hegger (SSH) coupling 1^. The local- 



ized spins correspond to the unpaired electrons of stable 
radical groups attached to the main chain, as realized in 
some substituted polyacetylenes 0, 0| . 

Thus, our model is presented in the Hamiltonian, 

H = -^t,(Ct,a+i^, + /i.c.) 

i i 

~'^iCl<7Ca ■ Sti + ClaC,2 ■ 8x2) , (1) 
ti = to + a{qi - qi+i) , (2) 



where Cj ^{Ci^s) creates (annihilates) an electron with 
spin s(=t or i) on site i. The nearest-neighbor transfer 
integral ti depends on the lattice deformation with the 
SSH coupling constant a. qi is the displacement of site 
i. U is the on-site Coulomb repulsion energy and rii^s = 
Cj gCi^s- K is the elastic constant of the lattice, whose 
kinetic energy is neglected. We define a dimensionless 
coupling constant A = /{IqK) and dimensionless bond- 
lengths = a{qi—qi+i) /Iq. The fourth term in the r.h.s. 
of Eq. (1) is the exchange coupling of two localized 1/2 
spins (Sti and 8x2) to the spins of tt electrons at site 
il and i2, respectively. C\ is defined as (C||,C||) and 
a = (fjj;, Oy^az) are the Pauli matrices. 

We exactly diagonalized the electronic part of the 
Hamiltonian in Eq. (1) by the Lanczos algorithm. Then 
the lattice is treated classically and optimized by means 
of the Hellmann-Feynman force equilibrium condition 
d (H) /dAi = {dH/dAj) = at zero temperature, re- 
sulting in 



A,^ xlj2{cia+i,, + h.c.) 



(3) 



where < • • • > denotes the ground state expectation 
value. 
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The alignment between the locahzed spins and the tt 
spins depends on the strength of the interaction parame- 
ters. Typical results are shown below with the following 
parameters in the unit of to (on the order of eV): U — 5.0, 
I J| = 1.0 while A is allowed to vary from to 0.5. The 
sign of J is taken to be positive (ferromagnetic) for most 
of the discussions hereafter, unless the antiferromagnetic 
(AF) coupling J < is specified for comparison. The 
origin of positive J has been discussed on the basis of 
molecular orbital calculations • 

Moreover spin alignment also relies on the length of 
the chain and the location of sites to which the local- 
ized spins are coupled. We calculated spin correlation 
function, total spin, and spin gap of the ground state for 
various chain lengths and spin locations in order to shed 
light on the spin alignment of molecular systems with 
different structures. In our model, the open boundary 
condition is imposed on the 7r-conjugated system. Since 
our model has the electron-hole symmetry, electron dop- 
ing gives the same results as hole doping. Thus we show 
the results for hole doping in the present paper. 

First we consider the case of the two localized spins 
coupled to electron spins at both ends of the chain, i.e. 
il — 1 and i2 = N. The spin correlation and the lattice 
deformation are plotted in Fig. 1 for both the half-filled 
case {Ne — N) and the single hole doped case {N^ — 
N - 1) with iV = 10 and A = 0.25 (a typical value for 
polyacetylene) . 

Figure 1 (a) shows the correlations of the localized spin 
Sti with TT electron spins and St2 • In the case of the half- 
filled system, the localized spins Sti and St2 exhibit an 
antiferromagnetic correlation mediated by the chain with 
the total spin 5 = 0. It is straightforward to understand 
this antiferromagnetic spin alignment for the half-filled 
ground state in terms of the topological rule. Because of 
antiferromagnetic correlation within the chain, any cou- 
ple of electron spins tend to have antiparallel alignment 
if there are an even number of sites between them. This 
spin correlation in the tt conjugated system together with 
the ferromagnetic exchange coupling J gives rise to the 
antiferromagnetic alignment between the localized spins. 

However, when a hole is doped into the tt electron sys- 
tem, the correlation between the localized spins becomes 
positive implying a ferromagnetic alignment. The spin 
correlation pattern in the right half of the chain (i > 6) is 
reversed by hole doping leading to the parallel alignment 
of the localized spins. Furthermore, the doped system 
turns out to be a spin quartet {S — 3/2). The inset in 
Fig. 1 (a) illustrates the alteration of spin alignment from 
the spin singlet to the spin quartet by such hole doping. 
This type of alteration corresponds to the one observed 
in the aforementioned thianthrene derivative 0. In the 
case of AF coupling (J < 0), single hole doping results 
in similar ferromagnetic alignment between the localized 
spins but with the total spin S = 1/2. 

The calculated lattice deformation patterns in Fig. 



1 (b) show that bond alternation appears in the half-filled 
system, while in the case of hole doping, bond alterna- 
tion becomes weak especially around the middle of the 
lattice. 

The electron-lattice interaction may be a two-edged 
blade to the spin alignment. On one hand, it can re- 
duce the spin correlation in the tt electron system because 
singlet pairs of electrons are formed in the dimerized lat- 
tice. On the other hand, the dimerization opens a Peierls 
charge gap to stabilize the system against the thermal ex- 
citation of TT electrons which may alter the spin structure. 

We have examined the effect of electron-lattice cou- 
pling in the range of A = 0.0 - 0.5. The spin correlation 
and the gap energy for both the half-filled case and the 
doped case are plotted as functions of A in Fig. 2. The 
absolute value of spin correlation actually increases as A 
increases for the half-filled ground state. For the doped 
case, the spin correlation is barely affected by A. 

Figure 2(b) shows the A dependence of the gap energy, 
and an inset indicates the total spin of the first excited 
state as well as that of the ground state. The gap between 
the quartet ground state and the doublet excited state in 
the doped case is slightly enhanced by the electron-lattice 
coupling, while in the half-filled case, the gap between the 
singlet ground state and the triplet excited state gradu- 
ally decreases with increasing A. Therefore these results 
show that weak electron-lattice coupling is preferable in 
the present case to achieve stable spin alignment control 
by charge doping. 

As is well known, 'tt topology' is of great importance 
to understand the spin alignment of organic molecular 
system. In the present model, it includes two major fac- 
tors: the size of molecules and the positions of localized 
spins. First we proceed to the size effect on the spin 
alignment. Figure 3 shows the correlations between the 
localized spins, the total spins, and the gap energies as 
functions of N with the localized spins located at both 
ends of the chain. For the half-filled case, the spin cor- 
relation follows the topological rule perfectly. The sys- 
tems with even numbers of sites are spin singlet {S = 0) 
with negative correlation, while those with odd num- 
bers of sites are spin quartet (S = 3/2) with positive 
correlation. Furthermore, even the hole-doped systems 
demonstrate regular behavior with respect to the num- 
ber of sites, although the topological rule is generally 
only valid for the ground state of the half-filled system. 
Thus hole doping changes the total spin in the following 
way: S = 0^ S = 3/2, if N is even; 5 = 3/2 ^ 5 = 0, 
if N is odd, as shown in the inset of Fig. 3(b). It is 
noted that the strength of spin correlation does not de- 
cay with increasing the chain length, although the gap 
energy gradually decreases as shown in Fig. 3(c). In the 
case of AF coupling, the total spin changes by doping as: 
5 = ^ 5 = 1/2, if TV is even; 5 = 1/2 5 = 0, if iV 
is odd. 

Next we investigate the spin alignment in the case of 
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the localized spins attached to inner sites of the chain. 
We pick up symmetric positions of the localized spins as 
{il,i2) = (l,Ar),(2,7V - I),---, for iV = 10. The spin 
correlation and the total spin are plotted as functions of 
il in Fig. 4. For the half-filled case, the spin alignment is 
antiferromagnetic with S = 0, as expected by the topo- 
logical rule, since there are even numbers of sites between 
the coupling sites. 

However, the localized spins in the doped case show 
positive correlation except for the location (il,i2) — 
(5, 6). Furthermore, those positively correlated spins can 
be classified into two categories in terms of the total 
spin. The doped systems are spin quartet (5 ~ 3/2) for 
il = 1 or 3 but spin doublet (5 = 1/2) for il = 2 or 4. 
A schematic picture in the inset of Fig. 4(b) illustrates 
the spin alignment for these two situations. These results 
show that the topological rule cannot be applied to the 
doped cases. The molecular systems with the spin loca- 
tions of (il,?2) = (1,10) or (3,8) are more suitable for 
the purpose of spin alignment control since they exhibit 
large change in total spin by doping. We have found that 
the violation of the topological rule also in the case of AF 
coupling. 

In the present paper we have shown the calculated re- 
sults for relatively large U and J, but even for much 
smaller values of these parameters the main results re- 
main qualitatively similar except the site dependence 
shown in Fig. 4. A thorough and detailed analysis of 
calculated results including such parameter dependences 
will be reported elsewhere. 

In summary, the spin alignment control in tt- 
conjugated molecular magnets by charge doping has been 
studied in the theoretical model of a Peierls-Hubbard 
chain connected with two localized spins. It is demon- 
strated that one-hole doping changes the total spin of 
the system from S — Oto5 — 3/2, if the localized 
spins are attached to the ends of the chain with even N. 
The intramolecular spin alignment is not much affected 
by electron-lattice coupling but depends sensitively on 
the topological structure of the molecular system. It 
is worth emphasizing that the topological effect in the 
doped case is different from that in the half-filled case. 
Conjugated molecules with controlled spin alignment as 
discussed in the present paper have a possibility of re- 
placing transition-metal ions in 'molecule-based' mag- 
nets such as Prussian-blue analogs, which exhibit pho- 
toinduced bulk magnetism JJ., 20]. 
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FIG. 1: (a) Spin correlations < Sti ■ Si > and (b) lattice 
deformations At for = 10 and A = 0.25. The half-filled 
{Ne — N) and the single hole doped cases {Ne = N — 1) are 
denoted by the dashed and solid lines, respectively. Inset in 
(a) shows schematically the change of spin alignment by hole 
doping. 



FIG. 2: (a) Spin correlations < Sti • St2 > and (b) gap 
energies as functions of A for the half-filled {Ne = N; dashed 
line) and the doped (A^e = N — 1; solid line) cases with A'^ = 
10. Inset in (b): the total spins of the ground and the first 
excited states. 



FIG. 3: (a) Spin correlations < Sti • St2 >, (b) total spins, 
and (c) gap energies as functions of N for the half-filled (Ae — 
N; open circle) and the doped {Ne — N — 1; filled circle) 
cases with A = 0.25. Inset in (b) shows schematically the 
spin alignment alteration by hole doping. 
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FIG. 4: (a) Spin correlations < Sti • St2 > and (b) total 

spins as functions of il, with i2 = A'^+l — il, for the half-filled 
(A^e — N; open circle) and the doped (A^e — N—1; filled circle) 
cases with N = 10 and A = 0.25. Inset in (b): schematic 
illustration of spin alignment with total spin S = 3/2 or 1/2 
in the doped case. 
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